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A mathematical model of the waterjet propulsion 
plant for surface effect ships consisting of gas turbine, 
reduction gear, axial pump, variable flushed inlet, diffuser 
and fixed nozzle was developed. Control algorithms were 
also developed for the variable geometry flushed inlet and 
for the pump RPM, and were tested by digital simulation of 
the mathematical model. Pump cavitation was included in the 
model and simulation. 

The control scheme for the variable flushed inlet 
involves the maintainance of a diffusion flow pattern about 
the flushed inlet. The control scheme for the pump involves 
the maintainance of pump RPM as dictated by the throttle RPM 
setting. In both control schemes, the rate and hard limits, 
on the hydraulic actuator for the flushed inlet and on the 
fuel flow rate of the gas turbine, were considered. Results 
of real time tests of the developed algorithm are presented. 
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INTRODUCTION 

The inherent geometrical flexibility of Waterjet 
Propulsion system commands a significant interest for its 
use in Surface Effect Ship (S.E.S.). In particular where 
the power level of the propulsion train is high, waterjet 
system, despite its relatively lower propulsive efficiency 
compared to the controllable reversible pitch supercavi- 
tating propellers (Fig. 1) is considered to be a more 
reliable system to fit the sidewall configuration of the 
S.E.S. A great risk is involved when propeller shafts at 
power levels of 20,000 H.P. and above are turned around 
corners to drive propellers. 

Thus the current U.S.A. effort for the projected 2,000 
ton S.E.S. having a cruise speed of approximately 80 knots 
focuses attention to the waterjet system as the best choice 
for propulsion. As a result of the experience of the 100 
ton waterjet driven S.E.S., it was recognized that a vari- 
able inlet geometry was needed to maintain high pressure 
recovery at all speed levels, particularly at hump speed and 
at cruise (Fig. 2) . Fig. 3 shows the typical speed resis- 
tance of the S.E.S. One of the design requirements is that 
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at hvimp speed the waterjet system must be able to deliver 
a maximum thrust of at least 20% greater than hump total 
resistance for an acceptable acceleration performance. 

This implies that the water flow rate for the system at 
cruise speed and hump speed will be approximately the same. 
Since the ratio of cruise speed to hump speed is approxi- 
mately 3, the flush inlet area must have to be adjusted 
accordingly so that a diffusion flow pattern will be main- 
tained always from the region in front of the flush inlet 
to the pvimp inlet (Fig. 4) . If not, then pressure recovery 
will be seriously degraded and the pump may cavitate due to 
insufficient suction head pressure. 

At present, work is being conducted at M.I.T. to 
optimize the design of components of the waterjet system for 
the S.E.S. so maximum payload can be obtained for a given 
displacement range , cruise speed and hump speed . 

The purpose of this paper then is to model and simu- 
late the dynamic performance of these conponents and to 
identify and bring out problems that may arise in the overall 
control of the system and possible criteria for the selection 
of the components. The case study presented herein is that 
of a 2109 ton S.E.S. whose waterjet system consists of six 
LM2500 gas turbines, reduction gears, axial water pvimps. 



6 







y 







iXl 

lU 



F16URG 4- OIFFO^U?Kl FtOW PATtcH»J AT IKILET AT HOMP AK/D CROl^jG 



7 



variable inlet geometry flushed inlet and fixed geometry 
nozzle. The simulation program has been set up in such a 
way that thrust reverser characteristics can also be added 
to analyze the backing performance of the craft. 
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CHAPTER I 

WATERJET PROPULSION COMPONENT MATCHING 
I-l GENERAL : 

In order for the waterjet system to function effi- 
ciently, pump design point must match accordingly to the 
thrust requirements both at hump and at cruise speeds. This 
includes the determination of the following: design range 

of variable inlet, area design pump RPM, design pump flow 
rate, design nozzle area. Using the momentum equation the 
thrust delivered by the system is: 

Td = ^ (Vj cosot - V^) (1-1) 

where: T^ = thrust develop (Ibf) 

= water density (Ibm/ft^) 

Q = water flow rate (ft/sec) 
g = gravity (32.2 ft/sec) 

Vj = waterjet velocity (ft/sec) 

Vq = S.E.S. velocity (ft/sec) 
oC = depression angle of the nozzle 
For most design cosot is approximately equal to unity 
thus can be neglected in the preliminary calculation. The 
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thrust develop then is: 




The thrust required at hiomp speed is: 



(1-2) 






'"SH + ■'lH> 



(l-3a) 



where: T ,, 

RH 



thrust required at hump 
speed (Ibf) 




bare hull resistance at hump 
speed (Ibf) 



R = inlet appendage drag at hump 
speed (Ibf) 



The factor 1.2 is the required minimum acceleration 



margin. 



The thrust required at cruise speed is: 



+ R-r^ 
RC SC IC 



(l-3b) 



where: T 

RC 



thrust required at cruise 
speed (Ibf) 




bare hull resistance at cruise 
speed (Ibf) 



Rjq = inlet appendage drag at cruise 
speed (Ibf) 



The inlet appendage drag for equations (1-3) is; 






(l-4a) 
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^IC “ ^DC ^OC^C 



(l-4b) 



where: C_„ and C„ = inlet drag coefficients 

DC a. u J 

at hump and cruise 
speeds respectively 

and = S.E.S, velocity at hump 

and cruise speeds 
respectively 

Qjj and = water flow rate at hump 
and cruise speeds 
respectively 



Equating thrust delivered and thrust required and 
substituting equations (1-4) to equations (1-3) results in: 



1»2 R3h 

-f VOH [ ^ -(1«-6 Cdh)] 

^sc 

-f ''oc [ ^ - <i+°-5=dc>] 



(l-5a) 



(l-5b) 



The pump head required at hump and cruise speeds from 
Bernoulli equation are respectively: 



2 

^PH " 2g 
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_ ^oc 

" 2g 




+ 2gH^ 
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(l-6a) 



(l-6b) 
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where: 07 = overall efficiency of inlet 

and diffuser 

= height of the nozzle above the 
waterline (ft) 



The effective horsepower required (EHP = p HpQ) for 
hump and cruise speeds from equations (1-5) and (1-6) are: 



EHPr 



EHP 

C 



1 . 2%Voh 1 




2 2gHwL _ m 1 

Vqh 2 toaj 


(2) (550) 


r VjH 
L VoH 


- (1+0.6 Cp)] 



RcVoc 1 




oal 


(2) (550) 1 




1 



(l-7a) 



(l-7b) 



Note that equations (1-5), (1-6), and (1-7) are all 

expressed as a function of jet velocity ratio (Vj/VQ) . In 
the case of a fixed inlet area the jet velocity ratio and 
flow rate at hump and cruise speeds can be related to the 
nozzle jet area (A^) by: 



_ Qh _ Qc 

^j “ VjH “ Vjc 



( 1 - 8 ) 



Fig. 5 shows a typical propulsive efficiency 
(RVq/EHP) versus jet velocity ratio (V^/V^) for different 
■ values of loss coefficient Kp where total head loss 
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«L = '/2g. Maximum efficiency is observed to occur 

for most practical designs between V^/V^ of 1.5 to 2.0. 
Likewise equation (l-7b) can be differentiated with 
respect to (V^/Vq) to obtain minimum EHP . However, the 
result may not be compatible to the requirement of equations 
(l-5a) , (l-6a) and (l-7a) because of the fixed nozzle area 

constraint unless the total hump resistance 
lower than some critical value. For the case where the 
total hump speed resistance is greater than this critical 
value the approach in finding the optimiam cross sectional 
area of the nozzle and jet velocity ratios at hxjmp and 
cruise speeds is as follows: 

a. Choose some multiple integral number of existing 
prime movers with a given shaft horsepower 
(SHP) rating. 

b. Multiply the maximum intermittent SHP available 
by the reduction gear efficiency and pump 
efficiency to obtain the maximum intermittent 
total EHP available. 

c. Substitute in equation (l-7a) and solve the 
quadratic equation for 

values of (V. /V^ ) represent the maximxjm and 
jH OH ^ 

minimum values of (V.„/V^ ) for which equation 

3H UH 

(l-7a) will be satisfied. 
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d. Solve for using equations (l-5a) and (1-8) 

using values of within the range 

obtained in step c. 

e. Solve for check equation (l-7b) 

if the total sustained EHP is within acceptable 
limits . 

f. Iterate steps a to e until optimum number of 
engines. A., propulsive efficiency is obtained. 

For the case where the total hump speed resistance 

is less than the critical value as dictated by equations 

(l-7a) , (l-7b) and (1-8) the procedure becomes easier since 

9 (EHPc) 

~ 3 (Vjc/ ' Voc) only to be set to zero to solve for 

optimum (Vj(^/VQQ) and the rest follows. 



1-2 PRELIMINARY PUMP DESIGN: 



The procedure stated herein applies primarily to 

determine the rated pump RPM, pump diameter, annulus area 

and other related pump characteristics based on designed flow 

coefficient (0) and head coefficient ('/') , designed flow rate 

(Q) and designed pump head (Hp) . By definition: 

Q 



^ ~ U “ NITD 



60 



(1-9) 
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( 1 - 10 ) 



where: C = water axial velocity (ft/sec) 

U = blade tip velocity (ft/sec) 

Q = water flow rate (ft/sec) 

2 

= pump annulus area (ft ) 

N = pump RPM 

D = pump tip diameter (ft) 



The axial pump annulus area (a) is: 



Where: 



A (1 - (d/D)2) 

n 4 

d = pump hub diameter (ft) . 



( 1 - 11 ) 



The design flow rate from equations (1-5) and the 
design pump head from equation (1-6) can be substituted 
in equations (1-9) , (1-10) and (1-11) to solve for Aj^, N, D 

for a given value of hub to tip diameter ratio. 



1-3 REDUCTION GEAR RATIO : 

The purpose of the reduction gear is to match the 
point of best efficiency of the pump to the point of best 
efficiency of the prime mover. Engine performance maps are 
usually expressed as plots of break horse power versus RPM 
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for various levels of specific fuel consumption (Fig. 6 and 
Fig. 7). Performance maps are also expressed as plots of 
torque versus RPM for various levels of fuel consumption 
(Fig. 8 and Fig. 9). In either case the point of best 
efficiency corresponding to power level required in equations 
(1-7) plus additional power for frictional losses and pump 
efficiency can easily be pinpointed. Thus the reduction 
gear ratio is determined as the ratio of the prime mover 
RPM at this point of best efficiency and the designed pump 



RPM (N) . 
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CHAPTER II 

MODELING OF S.E.S. AND WATERJET PROPULSION PLANT 
II-l GENERAL : 

In order to simulate the interaction of the com- 
ponents of the waterjet propulsion plant with the motion 
of the S.E.S. , it is necessary to obtain the mathematical 
model of such an installation. The model would then allow 
the development of the appropriate control algorithm with- 
out actually employing the hardware involved. To obtain 
such a model a thermodynamic approach was chosen. Fig. 10 
shows the power flow block diagram of the waterjet pro- 
pulsion train. At this point it is reconmended for the 
reader to take a look on ref. (1) and ref. (2) for the 
bond graft technique which will be used frequently in this 
paper . 




Figure 10. Block diagram of waterjet propulsion train 
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where: = fuel flow rate demand (Ibs/hr) 

= actual fuel flow rate (Ibs/hr) 

= power turbine torque (ftlb) 

^ = power turbine speed (RPM) 

= pump torque (ftlbs) 

Qj = pump speed (RPM) 

P 

Hp = pump total head (ft) 

T = thrust (lbs) 

Vq = velocity (ft/sec) 

The above diagram is interesting from the point of 
view of thermodynamic accounting of energies involved. The 
light broken arrows represent signal power level and the 
darkened arrows represent high power level transmission. 

The product of the variables affixed to the arrows multiplied 
by some constant represents power. For simplicity, the 
energy storage and energy dissipator elements associated 
with each component and power transmission are not included 
in the detailed modeling of each component which follows. 
Notice that all the energy associated with the fuel will 
in the end be dissipated into the environment primarily due 
to the S.E.S. resistance and secondarily but not avoidably 
due to the "resistances" inherent with the inefficiency of 
each component and inefficiency of power transmission. Thus 
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from the thermodynamic point of view, for the model to be 
valid in the steady state the power input for each com- 
ponent must be greater than the power output in either 
direction. Likewise the first law and second law of 
thermodynamics must be satisfied during the transient 
conditions . 

II-2 FUEL : 

The purpose of fuel is to provide energy primarily 
in the form of heat to the driving engine which is a gas 
turbine in this case study. The energy associated with the 
fuel can be expressed as heating value per unit weight (h) . 
Fuel reaction with air mainly is as follows: 

T 

FP 

2 „ + 21 12 CO + 18 H^O + 6 + Q, *• 

D lo 2 2 2 1 

12 CO^ + 18 H^O + Q;l + (2-1) 

where: T^p = flush point temperature 

Q^, ^2 ~ evolve in the combustion 

process respectively 

The degree on which the above combustion process proceeds 
to the right is dependent on combustion temperature, fuel-air 
ratio and time on which the process is allowed to proceed in 
the combustion chamber of the gas turbine. 




ii 
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Thus the total heat evolved per unit fuel weight as 
function of time (h(t)) can be assumed in block diagram 
to be: 



W. 







q2/qi 




O) 


1 +XiS 






1 +^2^ 


n 





where: q^, = heating values per unit weight 

of fuel from equation (2-1) 

h(t) = total heating value of fuel per 
weight 

= fuel flow rate 

= the exponential time constants of 
the above combustion process. 

X and X 2 very small terms but significant when con- 

vection time from fuel spray area of the combustion chamber 
to the turbine inlet is considered (Fig. 16). 

II-3 FUEL PUMP : 

The purpose of the fuel pump is to drive the fuel 
from the fuel tank into the combustion chamber of the gas 
turbine which is at a higher pressure. For the case of a 
positive displacement fuel pump (Fig. 9) represents the 



power flow diagram. 
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Figure 11. Pump power flow diagram 



where: 

cO 



R = 
m 



H 



op 

W, = 



= 



= 

m = 



input torque into the pump 
pump speed 

pump mechanical frictional resistance 
total lead developed by pump 
weight fuel flow rate 

fluid resistance associated with fuel 
value and fuel pipings 

pump torque 

fuel inertance associated in the fuel pipings 

modulo relating to H or modulo relating 

to ^ 



TF = ideal transformer of mechanical power to 
fluid power 

1 -junction = point of same CO or the same W^. 
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Since and are not independent for positive 
displacement pvunps, 1^^^ can be referred to the fluid side 
(or can be referred to the mechanical side) . Similarly 
can be referred to the fluid side. For the case where: 



(a) mW^ = 



or H 

op 




(2-2a) 



(b) = h CO 



(c) = c + 



Vlf 



f If 2gC2 (/® Aq) 



then 



(2-2b) 

(2-2c) 

92-3) 



, 2^ If, 

f"' Im T”’ TT 



= (mr.^ - m^bW^ - W^IW^I 



where: b = mechanical frictional resistance 

coefficient 

c-^ = fluid frictional resistance coefficient 
in pipings 

= orifice area of the fuel valve 
C 2 = orifice coefficient of the valve 

II -4 GAS TURBINE: 



Fig. 12 shows a schematic diagram of a typical air- 
craft gas turbine derivative. It is composed of two major 
units - the gas generator and the power turbine . The gas 
generator is a twin spool engine, consisting of low pressure 
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compressor (L.P.C.), high pressure compressor (H.P.C.), 
combustion section, high pressure turbine (L.H.T.) and 
low pressure turbine. The power turbine (P.T.) is adapted 
to the gas generator by an annular diffusing transition 
duct. Thus the operating speed of the power turbine is 
independent of the gas generator speeds. The engine also 
comes with a fuel control package to ensure that starting 
and operating schedules of the gas turbines are met auto- 
matically without the gas generator going into surge and 
without subjecting the turbine inlet to excessive gas 
temperature . 

In order to coordinate the gas turbine package to the 
water jet propulsion system, a mathematical model is developed 
as a part of this paper. 

For a satisfactory model, it can be assumed that the 
main accumulator of energy occurs in the gas pipe volumes 
and in the rotating parts. Although some heat energy is 
stored in the metal parts, its accumulation and discharge is 
carried out sufficiently slower. Thus during sharp load 
changes the transient is mainly from the dynamic interaction 
of rotors, volume and power turbine speed regulators. 

When a sharp load torque decrease occurs, the power 
turbine speed regulator initial action is to decrease quickly 
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the fuel flow to low value open up discharge value to 
release air intensively from gas turbine limit to evacuate 
the large quantity of stored energy. After which the 
discharge valve is closed and fuel flow is increased to 
the desired RPM. 

When an increase load is desired on the power turbine, 
it is necessary to increase the gas generator speed and 
therefore the stored energy therein. However in this 
transient process the compressor speed lags behind the 
pressure and temperature increase due to fuel flow increase. 
Thus not only for the protection of turbine against ex- 
cessive temperature but also preventing compressor insurge, 
the rate of fuel flow increase is strictly limited in 
accordance with air flow increase. Surge is also excluded 
at low speeds by discharging air at the first group of stages 
or from compressor delivery valve. 

In order to see the basic principle underlying surge, 
consider a multistage axial compressor having the same flow 
coefficient 0 = C /U and air axial velocity C for all 

X X 

stages, at design point. At very low speed C has to increase 
greatly through the compressor since the outlet air density 
is now approximately equal to inlet air density. This means 
that the first stages are operating at a very high incidence 
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angle or positive stall value while the last stages will 
be at negative incidence angle or negative stall value or 
turbining value. Now, with sudden increase of fuel flow 
the compressor outlet pressure will also increase very 
quickly. Thus C also decreases shifting the positive stall 
region toward the last stages. Since the first stages are 
not contributing any pressure increase, back air flow 
rapidly occurs as soon as stall hits the last stage with an 
accompanying explosive sound. 

With the above qualitative background. Fig. 13 is 
shown as the bond graph model of the gas turbine unit. 

Three variables are needed to describe power flow of com- 
pressible gas: mass flow rate (m^) , total pressure (P^) , 

total temperature (T^) . The subscripts are evident from 
Fig. 12. The rest of the notations are: 

= torque 

CO = angular speed 

= fuel weight flow rate 

R = frictional resistance for shafts or orifices 
I = polar moment of inertia 

H = total head of fuel 
o 

V = effective air volume capacitance 



FUEL 
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From the diagram the following are assumed: 



• 

in , 


• 

= m ^ , 


P ^ 


= p 


T 


= T , m 


• 

= m 


P = P , 


al 


a2' 


o2 


o3' 


o2 


o3' a3 


a4' 


o4 o5' 


T ^ 

o4 


= ^05- 


m 

a6 


= ^7' 


^o7 


0 

CO 

0 

P4 

11 


= T 

o8' 


^a8 = ^a9' 



P _ = P T _ = T m -_ = m TT 

o9 olO o9 olO alO all 



Above assumption implies that pressure drops along ducts 
and diffuser sections of the gas turbines must be included 
in the overall efficiency of the adjacent compressors and 
turbines . 

From Newton ' s Law: 



^L-cTE^ - ('^lt LC “ 

d CO H 

dt “ ^ ^HT ” ^HC ~ 

^P dt “ ^"^PT ” ^ PP ~ ^P^P^ 
Assuming air follows ideal gas law: 



P V = m CRT 
o o 



then: 



V 



1 



££o2 

dt 





) <RT 



o2 



+ m (RT 
1 o2 



dini 

= ”>a2 - 



(2-4a) 



(2-4b) 



(2-4c) 



(2-5) 



(2-5a) 



(2-5b) 
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where: 



where: 



dP 



V 



o4 



2 dt 



(iTi . - m — m ) (RT + m 

4a ad a5 o4 ,2 o4 



(2-5c) 



o6 



ad 
dP 



V_ = (ra _ - m _) (RT + m (RT 

3 at ' a5 a6 o6 3 o6 



(2-5d) 

(2-5e) 



^a5 ^dP^od^ 



o5 ^og) ^ 3 



“rTT “ IT ~ ^ 

at a5 a6 



(2-5f) 



<3Pq7 

V, — r— = (ifi ^ - m „) C^T „ + m,RT „ 
4 dt ' a7 a8 o8 4 o8 



(2-5g) 



dm4 

~7I = A ^ - ni „ 
dt a7 a8 



(2-5h) 



dP 



V 



olO 



6 dt ” "'alO^ ®'^ol0 ^5®'^ol0 (2-5i) 



dm 

5 • • 

= ra _ - m ^ ^ 

dt alO 



(2-5j) 



Equations (2-4) and (2-5) are the equations governing 

the dynamics of the gas turbine in Fig. 12. 

From compressor and turbine maps P ./P , versus 

o2 ol 

]^)/(Aj^Pq,) for different values of = — and polytropic 
^ ■'•ol 



34 




U>J e 6 OF o?P 5 TA K» T 
fpLViRopid. eFFfc-ieMc< 



C£>fJ5TA*JT ^ftE'P UN^6S 
TH(?U6APP5 Qf 






pl6Ugfcl-f cPft“g.AT)U6 KAf C?F c<PHf 




o S to t9 2S 




Ain ^/ 






0P£“A2AT/W6 A-<AP of AA/AL TURSft^e5 



35 



efficiency are plotted ^Fig. 14 and Fig. 15). Thus 

can be found as functions ^o2' ^ol' "^ol and annulus 
area (Aj^) . 

From Euler's equation, the torque on the rotor of 
an ideal lossless compressor or turbine is 

m^r ( ACe) 

g 

In terms of power: 



= m,<i> r A C„ = ra^ Aho 
a © a 



Therefore: 



m , Aho 



(2-6) 



Equation (2-6) will be true for an adiabatic compressor 
having no storage elements. Thus for a compressor having 
a polytropic efficiency 



^ polyc 



( y - 1) 
■ ■■ ' y 

(n - 1) 
n 



(n - 1) 



n 



therefore: 



(n - 1) ^ R 

polyc 
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where 



From equation (2-6) 



n (^-1) 

Po2~n~ 

Pol 


^02 


[^02 

S'^oiLToi ■ 


1 = C„T , 

J p oiL 




R 


vife) 




6) 






R 



(n- 1 ) 



Y__ = m C T . 
'^LC a p ol 



p<32\ 



f— ) 

\\^oV 



- 1 



(2-7) 

J 



(2-8) 



Likewise for adiabatic turbines having no storage 
elements and with a given ^poly,j,’ 



(n- 1 ) 

n 



^ PoIyt 



(n- 1 ) 

or -^- 7 ; — -RfV ^ 

^ ^POlYn 



therefore 



-'^LT "’aS^p'^oe 



^''Zpoly^ 



(^S) '• -J 



(2-9) 



Equations similar to ( 2 - 8 ) and (2-9) can likewise be 
developed for the H.P.C. , H.P.T. and P.T., thus defining 
the differential equations (2-4). Since from equation (2-7) 
the total temperatures and from performance maps mass flow 
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rates can be found and the rate of temperature changes can 
be reconstructed from the solution of the differential 
equations, equations (2-5) are then defined if the increase 
in total temperature inside the combustion chamber due to 
fuel energy is known. 

Fig. 16 shows the schematic diagram of the combustion 
chamber . 




Figure 16. Schematic diagram of the combustion chamber 



From the first law of thermodynamics: 



m^C T,. + W.Jc T + h(t)l = 
a5 Pa o5 fiNtpf o J- 



( 2 - 10 ) 



CW^ T^+C (W -W )T +C(m ^-m )T 
p fout o6 p^g fiN fout o6 p^ a6 a5 o6 



'Pa 



'Pfl 



^fg 
h (t) 



specific heat of air 

specific heat of liquid fuel 

specific heat of burned fuel 

instantaneous fuel heating value as 
defined in equation (2-2) 



where: 
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Since and X * 2 equation (2-2) are very small 

quantities, it can be assumed that h(t) = where h^ 

is the heating value of fuel per unit weight, since com- 
bustion chambers are usually designed to carry out complete 
combustion way up before turbine inlet. Thus, equation (2-10) 
in terms of fuel air ration W-/m ^ reduces to: 





o5 



+ 




Wf 

''Pfg 



+ c 



Pa 



( 2 - 11 ) 



is a significant term since this is the highest temper- 
ature which the turbine inlet is subjected to. For control 

purposes, the measured should account for: 

- s 

1. transport delay , e , where s is the 

Laplacian operator and is the time delay for the 

burning fuel to reach the turbine inlet, and 

2. time delay due to the inertia of thermocouples. 

It is worthwhile to mention here that among the 

quick acting means of temperature protection, fuel-air 
ratio limiters are widely used which are especially effective 
at starting and sharp load increase. However, fuel-air 
ratio limiters have low accuracy since changes in specific 
heat of fuel or changes in fuel gas ambient temperature and 
pressure do not make setting corrections . Thus there is 
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still a need for very quick acting thermocouples since the 
time constant for heating of the sharp trailing edge of 
turbine blades is in the order of one to three seconds. 

For the case of the LM2500 gas turbine, which has 
a single spool gas generator, equations (2-5a) , (2-5b) , 

(2-5g) and (2-5h) are simply deleted in the mathematical 
model . 

II-5 REDUCTION GEAR ; 

As mentioned in Chapter I, reduction gear is the 
impedance matching component to assure that the power 
turbine rotor of the gas turbine and the pump rotor 
operate with minimum fuel consumption at a steady state 
on power levels most often used. Thus it can be represented 
in bond graph model as in Fig. 17. 



I 



power TOKBiNE^ 








PUMP 



Figure 17. Reduction gear bond graph model 



40 



where : ^ 

COip 




m 



power turbine torque 

angular speed of power turbine 

pump torque 

angular speed of pump 

modulo of the transformer TF or the 
reduction gear ratio 



I = combined polar moment of inertia of 
pump and power turbine referred to 
the pump shaft 



= shaft bearing frictional resistance 



= reduction gear frictional resistance 
associated with reduction gear 
efficiency 



Assuming the reduction gear efficiency and shaft 
frictional resistance to be of the form: 



= 0.75 + 0.22(1.0 

L RG 




(2-12a) 



R^= C 2 + W (2-12-b) 

and the constants and C 2 are to be determined from 

data, then from Newton's Law: 



dCO 

I dt^ “ ^^T^RG - ^f “ ^p (2-13a) 




'I 










% 
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For the case where the power turbine is developing 
negative torque at idle, fuel flow rate during transients 
equation (2-13a) should be: 

<3 OJp m 71.J, 

I “HF- = - r (2-13b) 

P 



II-6 PUMP : 

Pumps are energy converters transforming mechanical 
power to fluid power. In order to understand the basic 
mechanism of pumps, consider the velocity diagram repre- 
sented in Fig. 18 for axial and/or centrifugal pumps 
without guide vanes. The solid lines represent the design 
point velocity diagram while the dotted lines represent the 
off design velocity diagram. 

U = rotor blade velocity 

= relative velocity of entrant fluid with 
respect to rotor 

= absolute velocity of entrant fluid to rotor 
and also exit fluid from stator 

C 2 = absolute velocity of exit fluid from rotor 
or entrant to stator 

W 2 = relative velocity of exit fluid from rotor 
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Figure 




For a fixed geometry pump, referring to Fig. 18: 



ACe = W^( cos - sin ^ tan f ) = 0{/S ) (2-14a) 

~ W^(cos^'- sin ^'tan ) = 0 (| 9 ') (2-14b) 



where : 






2 2 
W'^ = U' + C'^ 



sin /3 = ’ • 

\Tu2 + 

sin Q'= ^ ^ - -I 



From the above equation it can be seen that the 
change in the tangential velocity component ( of the 

fluid at any operating point is a strong function of 
dynamic similarity "flow coefficient" angle jS and a weak 
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function of angle "6 since the deviation angle of the fluid 
at blade trailing edges is approximately constant. Since 
sin and coSjS are periodic functions, then pump and/or 
turbine characteristics must necessarily also be periodic 
functions of dynamic similarity "flow coefficient" angle . 
From Euler's equation, torque ) and pressure (P) 

are: 



T' = p Qr 

P = p r AC oj 
> e 



where r is the pump radius 
where a? is angular speed 



Substituting equation (2-14) results in: 



= ]^p r0 ( ) ] . Q = (Rg ( ) . Q 
P = [pr0(|g ) ] • = <Rg( |3 ) .CO 



(2-15a) 

(2-15b) 



Thus, unlike positive displacement pumps and turbines, which 
are transf ormerly coupled, centrifugal and axial pimps and 
turbines are gyratorly coupled by Rg ( ) . 

If appropriate power dissipation terms, 
added to equations (2-15) such that: 



where: 



U-d = 



(2-16) 



subscript p = actual pump torque and pressure 



IPfj = power dissipation 
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then for a physical pump we have: 



tp 




■ Q + ra = Rg 


• Q + R^ • CO 


(2- 


-17a) 


'’p 


= (R 

g 




* • Q 


(2- 


-17b) 



where: 'C’d ~ dissipation torque 

= dissipation pressure 



IS : 



The bond graph from representation of equation (17) 




Figure 19. Axial or centrifugal pump bond graph model 

where: I = polar moment inertia of rotors 

m ^ 

1 ^ = fluid inertance in pipings 
= bearing frictional resistance 
R£ = pipe frictional resistance 

Note that the bearing frictional torque R^ is not the 
same as the pump torsional resistance R^ . By referring to 
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the typical centrifugal pump performance curves of Fig. 20 





Figure 20. Typical centrifugal pump characteristics 

and taking note of the points (x) where the power flow 't'CO 
or PQ are zero across the gyrator (GY ) , it is evident 
that the dissipation torque ^ ) and the 

dissipation pressure = P^ (Q, CO ). It can also be shown 
that and are not independent and can be related by 
some function R^(Q,CO). There the bond graph of Fig. 19 
reduces to Fig. 21. 




Figure 21. Axial or centrifugal pump model 
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where: R, = R.(Q,£0 ) = two port pump resistance 

a 



d 

R^ 



bearing frictional resistance 
frictional resistance on pipings 



From Newton's law and the above diagram: 



dco 

dt 



'^P - '^f - V 



(2-18) 



d^ _ 



•f dt 



«RgCO - R^Q) - Pf - Pp 



(2-19) 



where : 



'^p 

= 

^f = 
P = 



pump input torque 
shaft bearing frictional torque 
piping resistance friction pressure 
pump pressure load 



For S.E.S. water jet propulsion, the axial pump 

characteristic used in the system study (Fig. 22) was taken 

2 

from ref. (6), where the head coefficient = gHp/U is 
expressed as a function of flow coefficient C /U. Here 

the pump pressure of equation (2-17b) is: 

= — or P = = = (R„C0 - R O (2-20) 



where: p = density of fluid (Ih/ft^) 



g = gravity constant (32.2 ft/sec ) 
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FIOUQE 22 ?e,fifORMA^JCe CMAftAcTefJ.V^TvC of A/!.|AU 
PUMP 
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The design <f> and ^ of the above pump are 0.108 and 
0.51, respectively. To obtain "X ^ of equation (2-17a) , a 

ir 

dissipation function coefficient equivalent to equation (2-16) 
is assumed to be of the form: 

(2-21) 

min 

where: IP^ = dimensionless power dissipation 

IP^ = minimum TP , (at design point) 

“min “ 

<P^ = flow coefficient at design point 

= some constant to be determined 

relating to the efficiency charac- 
teristic of the pump 



Combining equations (2-21) , (2-16) and (2-20) , then 

t'p and Pp of equations (2-17) in terms of pump dimension 
reduces to: 

t'p = {^<P+ IPd) ^ [^'^] ^ (2-22a) 

P ^ Ictt ^ d1 

P = J (2-22b) 

P g 

where: = pump torque (ft-lb) 

ir 

= dimensionless head coefficient 

(p = dimensionless flow coefficient 

tP^ = dimensionless power dissipation of 
equation (2-21) 

P = water density (64 Ib/ft^) 
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= pump annulus area (ft ) 

N = pump RPM 

D = pump blade tip to tip diameter (ft) 



II-7 INLETS : 

Flushed and semif lushed inlets shown in Fig. 23 
and Fig. 24, respectively are considered to be more 
practicable to use for S.E.S. than ram inlets. 

The geometric parameters affecting the design 
pressure recovery and drag characteristics of flushed and 
semif lushed inlets are; entrance width to depth ratio, 
ramp wall divergence and ramp angle. For fixed geometry 
inlets, performance at various operating points is best 
described by efficiency of pressure recovery as a 

function of velocity ratio where is the velocity 

of the craft and is the average axial velocity of the 
fluid across the inlet throat area . Some typical pressure 
recovery performances of flushed and semif lushed inlets 
(ref. (32)) are shown in Fig. 25 and Fig. 26. 

Similarly, the drag coefficient 

Vi 

C = D_ ^ D . 

helo^i ■ Jf_ v^Q 

g 2g ° 
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FI6URE 23 FLU6HEl> IWLEf 




FI6URE 24 6EMJFLU6HE.0 IMLET 



Ram REcoveRv ratio 



51 




F\6URfc FUJ6HED IMuer RAM RECOVcRh' RATIO FOR FA«All£I- WAUL% 

I «• MeA^^uReo after iwlet section 

il- MEA60REO AFTeR pJFFO-5>£R 6EC.T10NJ ^ 

(WIDTH TO OcPTH RATIO = 2 *, RAMP AW6LE * '5 J 




FI6URE 0 .(» 6E Ml FLUSHED IWUET RAM RECOVERV RATIO 
I' MEASURED AFTER IWLET SECTlOW 
2“ MEA60RE0 AFTER P1FF05ER ^FCTIOU 



52 




Pl(aUKE 2^ PLO^HSO IWLET PRAGt 



la 

0 



a 

id 

O 

u 



2 

a 

o 



2 

oc 

UJ 

?< 

m 




figure 2<S pR^a OF GEMlFUJSHeCP tULET 



53 



can also be expressed as a function of the ratio V^/V^ 
to describe the drag (D) characteristics of the inlet 
at various operating conditions. Some typical drag 
characteristics of flushed and semiflushed inlets (ref. (32)) 
are shown in Fig . 27 and Fig . 28 . 

From the above pressure recovery characteristics 
the bond graph representation of inlet is: 







Figure 29. Flushed or semiflushed inlet bond graph model 



where: = fluid inertance in the inlet section 

= total pressure of the fluid way 
ahead of the craft 

= equivalent resistance of the inlet 
plus the boundary layer losses at 
the sidewall of the S.E.S. in front 
of the inlet 

“ total pressure at the exit of the 
inlet or at the entrance of the 
diffuser section so that at the 
steady state: P^^ ~ ^i^ol 
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II-8 DIFFUSER; 



The purpose of the diffuser is to increase the 
static pressure at pump inlet by converting the kinetic 
energy to static head. For a given expansion angle and 
fixed inlet condition, the total pressure loss of a 
diffuser is proportional to the theoretical value of the 

total pressure loss for a sudden expansion of the same 

• ^ou t 

exit area to inlet area ratio (A_ = ) . This pro- 

K ^in 

portionality constant is called diffuser loss factor (K^) 
and is defined by: 



K. = 



(^Ojn 



- P 



(^R - 



oput^^l 



(2-23) 



'R 



where : 



'-'in 

°out 

^1 

^1 

Ar 



total pressure at entrance to diffuser 

total pressure at exit to diffuser 

dynamic pressure at entrance to 
diffuser 

2g 

area ratio 



Figure 30 shows the diffuser loss coefficient as 
a function of area ratio (Aj^) and diffusion angle (20) for 
various geometrical shapes of diffusers. The bond graph 
representation, then, of the diffuser is: 



1055 FACTOR 
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pioufze 



PIPFO^BP- factor. 
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- fox I = ^3 

Q ^ Q 



R 



Figure 31. Diffuser bond graph model 



so that at the steady state: 



'’o 3 - "02 = 



(2-24) 



II-9 NOZZLE : 

The purpose of the nozzle is to convert the excess 
static head of the total head by the pump into kinetic 
energy head so that thrust can be developed. A well-designed 
nozzle, unlike diffusers, has negligible loss so that 
efficiency of 96-99% can be easily attained. The bond 
graph representation of nozzle is: 
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R 

Figure 32 . Nozzle bond graph model 

so that at the steady state: 

= "?n - V <2-25) 

where P^. = total pressure at the pump exit 

'-’in 

P = atmospheric pressure 

3 . 

= nozzle efficiency 

II-IO S.E.S., MASS, RESISTANCE AND WATERJET PROPULSION 
THRUST : 

The S.E.S. bare hull resistance data used in this 
study are tabulated in Appendix I. Note that if the 
boundary layer losses due to S.E.S. sidewall frictional 
resistance are included in the conputation for the total 
pressure recovered at the exit of the flushed or semiflushed 
inlet, then by energy conservation, an equivalent resistance 
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must be subtracted from the bare hull resistance data above. 
Otherwise, the propulsion system will be penalized twice. 

The assumed added mass for the forward motion of the S.E.S. 
was 10% of the total displacement. 

From the conservation of momentum, the instantaneous 
thrust developed for a fixed reference frame is: 



^ ^ /iT [4 J) [t{t ■ dS)] (2-26) 

V s 

or in terms of component thrust in x and y direction: 



/O Lj. dQ 

'^x = IT ^ 

pL dQ pQ 

T = — — + (V. sin o( ) (2-29b) 

y g dt g 3 

where: L and L are the distances between inlet 

X y 

entrance and nozzle exit in x and y directions, respectively. 

At the steady state, the "power input" to develop 
thrust is: 



Q 

in 




and power output is: 







(2-28a) 



TV 
d o 



(2-28b) 




i 
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where: 



APQ - T,V > 0 
d o ~ 



Since 





d g 



V. 

_JL 

V„ 



- 1 



then it is evident that the power conversion of equations 
(2-8) IS gyratorly coupled with modulo m^ - ^ 1'^ ” ^ 

and the bond graph representation is: 



1 5.6.^. 
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Figure 33. Bond graph model of thrust develop in S.E.S. 
water jet system 



The equivalent dissipator resistance 



Rd = 



AP- AP - AP 

loss _ in eff 



1 




I 
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where 



is: 



AP 



eff 



in 



1 

2 



5 L^o 



+ ni^ * V 
o t o 



R. 



\ fl _N_Zi _ Yi _ l] / Q 

g rol 1^2 IVol Vo Vo 2j / ^ 



(2-29) 



The overall waterjet system bond graph model which 
combines all the different components is shown in Fig. 34, 
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CHAPTER III 

SIMULATION AND CONTROL OF S.E.S. WATER JET 
PROPULSION PLANT 



III-l GENERAL : 

For the purpose of controlling the power turbine 
speed, the implementation of the detailed model of the gas 
turbine was found to be unnecessary. This was primarily 
due to the higher frequency response of the combustion 
process and pressures in the storage volume of the gas 
turbine in equations (2-5) compared to the frequency response 
of power turbine speed. However, for a very low value of 
polar moment of inertia of the power turbine drive train as 
referred to the pump shaft, the poles and zeroes of 
equations (2-4) and (2-5) will be significant in the design 
of a stable high frequency controller. Likewise, if the 
primary concern of designing a controller is the avoidance 
of surge in the compressor and of overtemperature in the 
turbine inlet of the gas generator then the detailed model 
will be necessary. 

In the case study it was assumed that the LM2500 
gas turbine used in the simulation was adequately protected 
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from overtemperature and surge by its own control package . 
This gas generator protective control setup was simulated 
by putting a scheduled limiter on the rate on which fuel 
flow demand signal can be increased. This rate can be set 
to have lower value than the built in limiter in LM2500 . 
The purpose of this limiter is to allow time for the gas 
generator rotor to build up speed and in turn air mass 
flow rate in order to neutralize the increase in pressure 
for surge avoidance and the increase in temperature due 
to sudden fuel flow rate increase. 



noted that the instantaneous torque developed at the power 
turbine is a function of the instantaneous air mass flov/ 
rate and power turbine speed. Air mass flow rate in turn 
is dependent on the instantaneous temperature and pressure 
at the storage volume of the duct between power turbine and 
the last stage of the gas generator. Thus the transfer 
function of the power turbine torque from equations (2-4 ) , 
(2-5), (2-7) and (2-9) can be approximated by: 



As for the rest of the gas turbine dynamics, it was 




(3-1) 



where: - = instantaneous power turbine torque 

£j 



^ss ” steady state power turbine torque which 
is a function of fuel flow rate (W^) 







I 
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and power turbine speed (N ) 

£j 



T’gg = time constant for the gas generator 
'C y = time constant for the volume between 



power turbine and gas generator 



s = Laplacian operator 



Since and are in the order of 10 ^ second, 

CjrVj V 



which is very much lower than the time constant of the 
power turbine drive train speed, the first term of equation 
(3-1) can be neglected so that: 



The plot of as function of fuel flow rate W 

^ ss f 

and the power turbine RPM is shown in Fig. 8. 

III-2 DYNAMIC EQUATIONS OF THE WATERJET PROPULSION PLANT : 

With the above approximation, the water propulsion 
system was reduced to the following non-linear ordinary 
differential equations: 




GG ^ i~+~~T 



1 



(3-2) 



a. 



The S.E.S. Velocity Equation 



dt W 




(3-3a) 
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where: = S.E.S. velocity (ft/sec) 

W = S.E.S. mass plus 10% added mass (Ibm) 

2 

g = gravity (32.2 ft/sec ) 

T = water jet propulsion thrust (Ibf) 

Rg = S.E.S. bare hull resistance (Ibf) 

Rj = flushed inlet appendage drag (Ibf) 

2 

= flushed inlet throat area (ft ) 



The instantaneous thrust (T ) in the longitudinal 



direction of the craft from equation (2-27a) is; 



T 

X 



g dt 




pQ 

+ (V. coso( - V ) 

g 3 o 



(3-3a-l) 



where: p = water density (Ibm/ft ) 

= longitudinal distance between the flushed 
inlet entrance and nozzle area (ft) 

= effective cross sectional area in 
longitudinal direction (ft^) 

= waterjet velocity at the nozzle exit 
(ft/ sec) 

= depression angle of the nozzle 



The S.E.S. bare hull resistance (Rg) should be a 
function of S.E.S. velocity (V^ and S.E.S. operating draft 
or air cushion pressure. In the simulation, however, it 
was assumed that the S.E.S. was at designed cushion pressure 
above hump speed. Below hump speed a resistance (R^) 
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proportional to was assvuned, so that: 



R3= 167.72 V„2 



0<V^50.7 (3-3a-2) 

o 



R = 26.726 V - 6935.34 V + 712204.0 V>50.7 

so o o~ 



The flushed inlet drag R from equation (2-23) is: 

V, 



Rt = ■“ Cr, — Q ^ 



'I ■ 2 D 



g ° ' Vq 



(3-3a-3) 



For the purpose of controller design the flush inlet 
drag coefficient was assumed to have a sharper point of 
minimum value at Vj^/VQ=0.8. The minimum value of was 
also assumed to be equal to 0.15 consistent with the published 
experimental data, so that: 

Cd = 1.0 

C(3 = 1.0 - (0.8/0.85) (V^/V^) 1.0<V^/V^0.8 

C, = 0.2 V-/V - 0.01 V./V^>0.8 

d 1' o 1 o 



b . The Pump Speed Equation 



i? = f^[t3(r^,RGR) - -!rp(N) - rp(N,Q)] o-sw 

where: N = pump RPM 

= polar moment of inertia of power 

turbine, reduction gear and pump ^ 

rotor referred to pump shaft (ft-lb/sec ) 
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= effective engine torque at pxomp shaft 
(ft-lbf) 

7T p = shaft bearing friction torque (ft-lbf) 
= pump torque (ft-lbf) 



The effective engine torque 









equation (2-13) is: 



n = '^E 7rg 



Tc. = (RGR) 



'^RG 



XT >0 

E— 






where: RGR 

<^RG 

^^RG 



(3-3b-l) 



= reduction gear ratio = 5.55 

= reduction gear efficiency 

-Nt /895 

= 0.75 + 0.22 (1 - e ^ ), where 

Nc. = power turbine speed 



The shaft bearing torque ( Z^„) is assumed to be: 

F 





= 1000 + 
F 


15N 


(ft-lbf) 


(3-3b-2) 


The pump 


torque ( Xp) 


from 


equation (2-22a) 


is: 


r 

p 


= {^4 + 


2g 


-f 2 

[;•] 


(3-3b-3) 



where 

_Q 
An 

= N1?D 

~EU~ 
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A curve is fitted to the pump characteristics of 
Fig. 22 for the head coefficient ( 'Z' ) as a function of 
flow coefficient (^) so that: 



-83. - 1.916^ + 0.85 
\j/= 220(f>^ - 25.5(p + 1.308 

-68.217 + 10.173^ + 0.2069 

-318.45^^ + 72.713^ - 3.6949 



^:< 0.04 
0.0 < 0 ^ 0.065 

0.065< ^<0.130 
0.130^^0.152 



The dimensionless power dissipation tP^ of equation 
(3-3b-3) was normalized to a pump having an efficiency of 
90% at design point so that equation (2-21) yields; 

= 0.00623 + 12.4 ( 5^ - 0.108) 

c . Waterjet Flow Rate Equation per Pump 



dt 



D 



Hp(N,Q) 



H (V^,A.) 
SV o i' 



Hj (Q) 




(3-3c) 



3 

where: Q = water flow rate (ft /sec) 

_ 3 

P = water density (Ibm/ft ) 

2 5 

= fluid inertance in pipe (lb-sec /ft ) 

Hgv = pump net positive suction head (ft) 

= waterjet kinetic head at the nozzle (ft) 

= head generated due to S.E.S. 
acceleration (ft) 
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• *, 

The pvimp head from equation (2 -2 2b) is: 

r N 1 2 

Hp = g[^^Dj (3-3C-1) 

The net positive suction head from equation (2-25) 

and from flushed inlet efficiency characteristics is: 



H 

SV 


2 

= ? ^ ^i"?D - Hwl (3-3C-2) 


where: ~ 


flushed inlet efficiency 




diffuser efficiency 




pump height above waterline 


7i = 




= 





The water jet kinetic head at the nozzle (Hj) is: 





2 

""j 

“j = 2g(^N (3-3C-3) 


where: cn 

7 N 


= nozzle efficiency “SJ 1.0 


''j 


= Q/C^Ajjoz ~ waterjet velocity (ft/sec) 


Cc 


= nozzle coefficient of area con- 
tractions 1.0 


^OZ 


= cross sectional area of nozzle 
throat (ft^) 
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The head generated through S.E.S. acceleration (H^) 



is: 






Lxf^g ^ ^ 

= g [dt - dt__ 



(3-3C-4) 



where: L = distance between entrance to flushed 

X 

inlet and nozzle exit in the V 

. . / o 

direction (ft) 



d . Power Turbine Torcfue Equation 

= 0.5 \Z ) - 't' 1 

dt I ss f E eJ 



(3-3d) 



where: = instantaneous power turbine torque 

(ft-lbf) 

^gg = steady state power turbine torque 

as found in the engine performance 
map 



N 



g = power turbine RPM 



The factor 0.5 is the ( term in equation (3-2). 



e . Fuel Flow Rate Eguation 
dW^ 

= 2.0 (W^^ - W^) (3-3e) 

= instantaneous fuel flow rate (Ibm/hr) 

= fuel flow rate demand signal (Ibm/hr) 



where : 




J 
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Equation (3-3e) assumes that LM2500 control package 
has a built-in servomechanism for the fuel pump with a 
time constant of 0.5 second and a constant gain of 1.0. The 
rate limit and hard limit were incorporated in the fuel flow 
rate demand signal so that 



"fd ^ "f d. 



max 



W 



fd, 



min 



4L < W 

- fd - fd 



max 



The minimum fuel flow rate demand signal was 

fd • 
min 

set to idling fuel flow rate of 800 Ibm/hr and the maximum 

fuel flow rate demand signal was set to 10,000 Ibm/hr, 

^°max 

The maximum rate of change of the fuel demand signal 
was set to six seconds from to W_, 



^max 



min 



max 



f . The Flushed Inlet Area Servomechanism Equation 



dA^ 

dt 



=0.5 (Aid - A^) 



(3-3f) 



where: A. = instantaneous flushed inlet throat 

area (ft ) 

, 2. 

A = flushed inlet throat area demand (ft ) 
id 

Equation (3-3f) assumes a hydraulic servoactuator for 
the flushed inlet throat area with a time constant of two 
seconds and a constant gain of 1.0. In the simulation, the 
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2 

hydraulic actuator was set to have a hard limit of 7.0 ft 

2 2 , 
maximum, 1.5 ft minimum, and a rate limit of 0.5 ft /sec 

so that 



A. < 0.5 
1 

1.5 < A£ < 7 .0 
III-3 PUMP CAVITATION: 



To simulate pump cavitation during transient 
condition at various pump operating points, a dimensionless 
suction head criteron was assumed to be: 






SVC 



(<^) 

' SVC ' 



where: „ = critical net positive pump suction 

head (ft) 

U = pump blade tip velocity (ft/sec) 



A curve is fitted into the pump cavitation data 
(Fig. 22) so that is defined as a function <f) . The 

pump is assumed to cavitate if the actual dimensionless 



suction head ^ 5 ^ 



SVA 



U 



sv fell below 



SVC 



The degree 



of degradation in pump performance was assumed to be of 
the form: 

2 



\JJ 1 - r r ^SVC~ ^ SVaI 

^C'T -^i[ J 



V^SVA ^ ^SVC 
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where: '^q ~ dimensionless head with cavitation 

dimensionless head without 
cavitation 

C-^ = some constant to fit the pump 
cavitation data 



Substituting equation (3-5) to equation (3-3c-l) , 
the pump head generated with cavitation is: 



H 









SVC ^ SVA 






SVC 



[eo i^svA<fsvc 



(3-6) 



The dimensionless power dissipated when there is 

cavitation [P, is also assumed to be of the form: 
dc 



tP = p 
dc a 



1.0 + C, 



[ 



^ SVC 'i' SVA 



f 



SVC 



] 



^SVA <*/^VC 



(3-7) 



where: (P^ = dimensionless power dissipated without 

cavitation as defined in equation (2-21) 

C 2 = some constant to be determined from 
pump data 



Substituting 'p r, and P, in equation (3-3b-3) , pump 

dc 

torque with cavitation is: 



^ ‘fsvA<V'svc 









‘f’ 





Nil. 



’-=1 



!1 
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III-4 LINEARIZED DYNAMIC EQUATION OF THE WATERJET PROPULSION 

The linear expansion of equation (3-3a) about a 
steady state operating point yields: 



V = 






3T 

^ <Voo>'" 






(Q_Q ) S. (V-V ) 

9 Q ^o^ -av o oo^ 



0Ri 3R. 1 

bV (%^o<i “ 3q “ 0Aj^ ^^i“'^io^J 



The subscript ^ represents the steady state operating 
point values of the state variables and the partial deri- 
vatives are evaluated at the corresponding steady state 
operating point values. Similar expressions can be derived 
for the rest of the equations of (3-3) . At steady state: 



Voo = = Qo = = *io ' ° 



Defining : 



V-V = V 
o oo o* 



and 



Vo-Voo = Vo = V 



O* 



Q-Qo = Q* 



and 



Q-Qo = Q = Q_ 



etc. yields the linear form of equations (3-3). 
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^33 



If L 3 q-* 



A 



36 " 



A 



42 



A 



44 





0.5 RGR 






SS 



dNr 



A 



55 



2.0 




- 2.0 



A 



66 



0.5 




-0.5 



B 



51 



2.0 



B 



62 



= 0.5 




2 .0 



0.5 



The block diagram of equation (3-8) is shown as 
solid lines in Fig. 35. The broken lines represent the 
controller. The cuts in the diagram at points x , y , 
and z will be discussed in the succeeding sections. 
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III-5 FLUSHED INLET THROAT AREA CONTROL ; 

For the performance of semiflushed or flushed inlet 
to be optimum 1) pressure recovery must be maximum and 
2) drag must be minimum. Fortunately, these two criteria 
were not contradictory, as can be seen in the performance 
curves of Fig. 25, 26, 27, 28 of Chapter II. At the inlet 
to ship velocity ratio of ~ 0.8, the above two con- 

ditions are satisfied. A control algorithm then can be 
set up in such a way that equivalently 

O/Aj^V^ = 0.8) at steady state and at least close to 0.8 
during transients. A preliminary proportional and integral 
control for the flushed inlet actuator was set up in the 
form of: 




(3-9a) 




(3-9b) 



so that: 




(3-9c) 



Substituting (3-9) to (3-3f) yields: 



A. 

1 



-0.5Aj_ + 0.5A 




- 0.8) (3-10) 
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Taking the linear form of equations (3-9b) and 
(3-10) and including them in equations (3-8) yields; 



V* 




1 — 1 
1 — 1 


0 


^3 ° 


0 


^16 


0 








0 


N* 




0 


^22 


^23 ^24 


0 


0 


0 




N* 




0 


Q* 

* 




^1 


^32 


A33 0 


0 


0 


0 




Q* 




0 


T'e* 


= 


0 


^42 


0 

44 


^45 


0 


0 




E* 


+ 


0 


Wf* 




0 


0 


0 0 


^5 


0 


0 




^f* 




®51 


• 




\l 


0 


\3 ° 


0 


^66 " 


‘67 




^i* 




0 






^71 


0 


A 73 0 


0 


^76 


0 








0 


where: 










1 












— 






^UiVo 


1 






Q 








- U._) 


a Vo 






^A^A 


. V ^ 














1 












^63 








j 






1 1 






- U.D 


3Q 




• 1 


PAIa. 


vj 



] 



(3-11) 




A = 
66 



Ag6+0 • 5Kp;^— 








= 0.5 




-K 



lA 



[ 



A- V 2 
o 



] 
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^3 = 



lA 



[^] 

lo 



A = -K,, 
76 lA 



[rfr] 



1 o 



This PI control is shown in Fig. 35, in which; 



Bi = AiVo^ 



Bo = A. 

2 1 o 



Bo = — S' 



3 V Vo 



To obtain the steady state values of V , N , Q , 

oo o o 

/ W, and A. an open loop response to different 
"o ^o 

levels of using conservative values of K and K 

fa lA PA 

were obtained. Some typical responses are shown in Fig. 36. 

To obtain an optimum K and K , an open loop 

lA PA ^ 

transfer function by cutting the loop at point x 

in Fig. 35 was obtained using MBODE program. Such open 
loop transfer functions can be represented in matrix form: 
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V* 




^1 


0 


^13 


0 


0 


0 


N* 




0 


^2 


^23 


^4 


0 


0 


Q* 




^31 


A 

32 


^3 


0 


0 


0 


E* 


= 


0 


^42 


0 


44 


^45 


0 


W 




0 


0 


0 


0 


A 

55 


0 






^61 


0 


^63 


0 


0 


^66 


• 

A. , 




^71 


0 


A 

73 


0 


0 


^76 


[v] 


= 


[o 


0 


0 


0 


0 


1 



0 

0 

0 



0 



“] 





V 

* 




’ l' 




N 

* 




0 




Q* 




0 






+ 


0 




Wf 




0 




A. 

1* 




0 


J 






0 



[’•] 



V* 

N* 

Q* 

t 

w. 



E* 



1 * 



id. 



where of equation (3-11) is set to zero. 

The zero location can be adjusted by the ratio of 
Kj^/Kp^ depending on the location of poles and zeroes of 
the hydraulic servoactuator for the flushed inlet at various 
operating points. The hard limit and rate limit of the 
actuatoi; on the other hand, sets a constraint on the 
maximum value of Kj^ that can be used, otherwise excessive 
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overshoot occurs. For the case study presented, the hard 
limit and rate limit on the actuator was overcome by 
setting = 0.5 and Kp = 5.0 and putting a hard limit 
on the integrator portion of the controller equal to the 
hard limit of the actuator, so that the maximum and 

minimum values of A. . were 7.0 and 1.5, respectively. 

I 

Thus for the control of the flushed inlet geometry 
some means to measure Q, V, must be provided in order 
to keep track of V^/V^. The time delays associated with 
these measurements must also be included for the proper 
location of the zero of the PI controller. 

III-6 OPEN LOOP SIMULATION TO FUEL DEMAND SIGNAL ; 

The fuel demand scheduled limiter used in open loop 
simulation was set to have a maximum rate of change from 
idling to maximum in six seconds and a hard limit 

of idling = 800 Ibm/hr and maximum = 10,000 Ibm/hr. 



term in the thrust equation (3-3a-l) were neglected since 
they were negligibly small compared to the other terms. 
Since S.E.S. under study uses six engines, the thrust 
equation (3-3a-l) and the inlet appendage drag equation 
(3-3a-3) were multiplied by six for the simulation of the 



The term in the flow equation (3-3c) and the 



o 
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straight line motion of the craft. Computer runs were 

made using a stepped fuel demand signal from each 

engine idling fuel flow rate of 800 Ibm/hr to various 

levels of up to maximum of 10,000 Ibm/hr per 

engine. The stepped W signal was filtered by the maxi- 

fd 

mum rate change of the Scheduled W limiter with the 
setting discussed above. For a realistic open loop simu- 
lation, the flushed inlet area control loop discussed in 
the preceeding section was closed. Pump cavitation was 

observed at the stepped W of 5,000 Ibm/hr/engine and 

f d 

above. Some responses of pump RPM, pump flow rate (ft /sec) , 
S.E.S. velocity (ft/sec) and flushed inlet area (ft ) to 
the stepped made are plotted in Fig. 36. 

Fig. 36 responses, however, were obtained using a 
higher value of fluid inertance (I^) in the pipings than 
the computed value for a faster YSYS solution (i.e., using 
higher value of time step) . Comparison of the response with 
the actual value of I (i.e., using very small value of 
time step) yields a similar result. The higher value of 
l£ used served as a filter to the oscillatory solution of 
YSYS to pump flow rate during transient and especially in 
the cavitation regime when the time step is greater than 



0.1 second. 



350 
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III-7 PUMP RPM CONTROL: 



Since pump shaft RPM appears to be the easiest to 
measure accurately and appears to be the most effective 
indicator of power turbine stress overload, it was decided 
to base the fuel demand signal on the error (Ejj) between 
the throttle RPM setting (Nj^) and the actual pump shaft 
RPM (N) . 

By inspection of Fig. 35, the open loop transfer 

function yields a type O system. Such open loop 

transfer function, represented in matrix form of equation 

(3-11) with the output matrix [c] = (]010000 0][ , was 

verified using MBODE program and the result yields no pole 

at the origin in the s plane. Thus, to reduce the error 

(Eu* = - N^) to zero at steady state, an integral 

control is needed . A proportional control is also needed 

to bring close to faster for large error (Ejj) . In 

order to obtain the gains for above PI control, 

MBODE program was used to locate the open loop poles and 

zeroes of for various levels of i.e., 

^ rci* rd 

800 Ibm/hr to 10,000 Ibm/hr. 



In the MBODE solution, however, of equation (3-11) 



the A 23 term was so big compared to the ’^ 2 '!' ^ 44 ' ^55 



and 



'66 



terms that a big numerical error resulted. It 
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was therefore necessary to reduce the order of the system 
to six by setting: 



Q* = 0 



V + A 
31 * 32 



N + A 
* 33 



Q 

* 



This is equivalent to setting the fluid inertance to zero 
or shifting the pole due A^^ to - oo in the s-plane. Such 
a procedure will yield a solution very close to the true 
solution since the pole due to being so far away to 

the left of the other system poles and zeroes, is effect- 
ively at - oo . Solving for Q yields: 

* 



Q* 





(3-13) 



Substituting equation (3-13) into V^, 
of equation (3-11) results in: 




and A 



id. 




^1 


^22 


0 


0 


a 

15 


^21 


^22 


^23 


0 


0 


0 


^32 


^33 


^34 


0 


0 


0 


0 


^44 


0 


^51 


^52 


0 


0 


^5 


a 

61 


^62 


0 


0 


^65 



0 




V 

* 




‘ o' 


0 




N 

* 




0 


0 








0 


0 




"f. 






^56 




A. 

1* 




0 


0 




id^ 




0 
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where : 



^11 = -^11 



^13 ^31 
^33 



22 



^13 ^32 
^33 



^15 “ \e 



*21 



^3 ^1 

^33 



^22 ^22 A 



^23 ^32 



33 



^23 ^24 



332 - ^42 



333 - A 



44 



^34 “ ■^45 



%1 



= A 



^63 ^31 



61 A 



33 



^2 “ A 



^63 ^32 



33 
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With the above PI control, since W = 2.0(W^,-W ) 

f fa f 



K 



="'5 "fd = 



IN 



g ) (N -N) , then the closed loop 
R 



transfer function N /N is represented by: 

* R* 



V 


' 




a, .. 


0 


0 


ai r- 


0 


o" 








o' 


* 




11 


12 






15 








* 






N 

* 




^21 


^22 


^23 


0 


0 


0 . 


0 




N 

* 




0 






0 


^32 


^33 


^34 


0 


0 


0 








0 







0 


a . « 


0 


a 


0 


0 


a , 






H- 


b. 


f * 






42 




44 






47 




f* 




4 


^i* 




^51 


^52 


0 


0 


^55 


^56 


0 




A. 

1 * 




0 


• 




^61 


^62 


0 


0 


®65 


0 


0 




A. . 

id^ 




0 






0 




0 


0 


0 


0 


0 
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. j 


with the output matrix 


W 


= [o 


1 0 


0 0 


0 


o] 





(3-15) 



N 



R 



where : 



a ——2.0 Kjju 
42 PN 



a =+2.0 
47 



^7 2 *^IN 



b = +2.0 K 
4 PN 



^ =-^^IN 
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a = A 

55 -^66 





A 7i 


= A 




71 




33 


^73 


^32 



Vn = — A N* 

62 -^33 



"65 = ^6 



i> = B 
4 "l5 



An MBODE result for the value of matrix A of 
equation (3-14) at = 5,000 Ibm/hr with output matrix 

£c]] =[ 010000 ] indicates that three of the zeroes 
cancel out three of the poles of the open loop transfer 
N^(s)/W^^ (s) and is plotted in Fig. 37. 




Figure 37 . Open loop poles and zeroes of N VWj. , 

* fd* 
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Thus in the PI controller of the form: 



^IN ^PN 

+ = 

PN s s 



K 



IN 






+ s 



N 



Ki^/Kp^ can be adjusted so that the zero of the controller 
is in the vicinity of -0.234 and the gain also adjusted 

so that closed loop pole is within the 45° line indicated 
in Fig. 38 for a fast response v/ith less than 5 % overshoot. 



/ 







PN 



Figure 38 . Root locus plot for 



0.234 




V- 
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III-8 ANTI WINDUP LOGIC FOR PUMP RPM CONTROL ; 

In order to compensate for 1) the excessive over- 
shoot and oscillation during transient caused by the hard 
limit and rate limit on W^^ signal and 2) the shifting of 
the two rightmost uncancelled poles of equation (3-14) 
towards the origin at low power level, the K and K 
were adjusted accordingly in the form as shown in Fig. 39. 

The idea behind such a scheduled and Kpj^ is to 

set the gains of the PI control for small error signals to 
be equal to the optimum values as found in the linear 
analysis of section III-7, so that: 






max 



K = f (N ) 

^max 2 R 



The f ^ ^2 plots of the optimum 



K 



PN 



and K as function of N which will yield a fast stable 

response with less than 5 % overshoot about an operating 

point N^ for small error disturbance. For large error 

disturbance, however, the value of K is too large to 

^max 

be compatible with the rate and hard limit of the scheduled 
fuel demand limiter. Large overshoot occurs because by the 
time the integral part of the fuel demand signal ) has 



96 



KfhJ 




Kin 
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reached the correct value of W for a steady state N , 

f R 

the error signal is still lagging too far behind. This 
is due to the fact that the correcting signal acting 

on the fuel pump servo is considerably less than the 
combined integral and proportional fuel demand signal 
(W_ , + T ) due to the filtering action of the scheduled 

fdj fdp 

limiter. Since for large error the proportional part 

of the fuel demand signal ) by itself is large enough 

to be filtered by the limiter, therefore can be 

set at a lower value and still get the same speed of 

response in pump RPM (N) . Thus K_ can be set so that 

■‘•min 

by the time the error signal goes to zero, the integral 

part of fuel demand signal (W^j ) is equal to the actual 

fuel flow rate (W ) for reference pump RPM (N„) . Over- 

r K 

shoot through this setup is eliminated . 

However, for IC = constant = K , it was observed 
^ ^max 

that some secondary oscillation occurs during the accelera- 
tion transient (Fig. 33b). This is due to the hard limit 
on the signal for large error. When the error signal 

v/as decreased to a value such that , was no longer hard 

fd 

limited, the W signal decreased rapidly, which caused 
f d 

the actual fuel flow rate to decrease slightly. To 
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overcome this effect, a negative sloping function of Kpj^ 

with respect to absolute value of error was set so as to 

maintain always an increasing fuel flow rate (W^) during 

acceleration transient. K was also set so that for 

^min 

large errors will still be hard limited - 

III-9 CONTROL ALGORITHMS TEST RESULT : 

The following hard maneuvers were simulated to test 
the control algorithms for pump RPM and for the flushed 
inlet throat area: 

(a) Acceleration from 109 to 635 RPM Throttle Setting 

This simulation was conducted by putting a step change 
in RPM throttle setting. The 109 RPM corresponds to idling 
power setting and the 635 RPM corresponds to 95% of the 
maximum intermittent power rating of the engine. The results 
are plotted in real time (seconds) in Fig. 40. The labels 
and scales are: 

V = S.E.S. velocity (ft/sec) 

N = pump RPM 

3 

Q = water flow rate per pump (ft' /sec) 

Number of pumps = 6 

WF = actual fuel flow rate per engine (lbm/hr)X100 .0 



Number of engines = 6 
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AI = actual flushed inlet throat area per inlet 
(ft2)X50.0 

Number of inlets = 6 

TRST = total thrust of the six (6) waterjet system (Ibf) 

TE = engine torque (ft-lbf) 

TP = pump torque (ft-lbf) 

Notice that the gas turbine is dually protected from 
surge, overtemperature and transient mechanical overstresses 
by 1) the fuel demand rate limiter, and 2) by the PI con- 
troller gains scheduler of Fig. 39. Thus for large errors 
the fuel demand signal received by the engine will at most 
be a six second ramp signal. This is verified by the WF 
plots in Fig. 40. 

The response in pump RPM shows that the controller 
setup of Fig. 39 is fast and stable. In addition, it has 
no overshoot. There is also no secondary oscillation due 
to the rate and hard limits of the scheduler during 

transient as this was taken care of by the negative sloping 
function of Kpj^ with respect to the absolute value of error. 
Thus as the error decreases the fuel pump still receives an 
increasing signal. 

The flushed inlet throat area (A ) increased in the 

i 

initial transient to accommodate more flow at low S.E.S. 
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speed to generate the acceleration force. It finally 
decreases as the S.E.S. builds up speed, maintaining 
close to 0.8 value during transient and exactly 0.8 at steady 
state. Thus, setting a hard limit alone on the integral 
part of the PI controller for the flushed inlet is enough 
to overcome the rate and hard limit of the flushed inlet 
actuator . 

The pump torque ( 't ) and engine torque { 'C ) are 

P E 

plotted to check mechanical stress overloads during 
transients . 

(b) Deceleration from 635 to 150 RPM Throttle Setting 

This simulation was also conducted by putting a step 
change in RPM setting. In order to protect the gas turbine 
gas generator turbine inlet from sudden change in tempera- 
ture the fuel demand rate limiter for decreasing fuel flow 
rate was set arbitrarily to six second ramp. The PI con- 
troller gains scheduler of Fig. 39 likewise contributes to 
this protection as an added precaution. The results are 
plotted in Fig. 41. The labels and scales are the same as 
those of the acceleration test. 

From the plots it can be seen that the response of 
pump RPM was slower than in acceleration and has slight 
overshoot. This was expected since the damping term in 
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pump RPM varies approximately to the square of pump RPM. 

Such an overshoot, however, is within acceptable limits 
and in addition it occurs at approximately 10% of the 
maximum intermittent power rating of the engine. 

Thus, it is evident from the pump RPM responses 
that the PI controller gain of Fig. 39 works optimally 
throughout the whole operating range of the waterjet system, 
that is, 1) fast response, 2) less than 5% overshoot, 

3) stable and 4) no steady state error. 

From the plot of it can be seen that the actual 
fuel flow is rate and hard limited even without the actuator 
itself being rate and hard limited. 

The plot of the flushed inlet throat area (A^) showed 
an initial decrease to its minimum hard limit for several 
seconds in an attempt to bring equal to 0.8. This did 

not present a problem on the integral part of the PI con- 
troller for the flushed inlet since the integrator was also 
set to have a minimum hard limit equal to the minimum hard 
limit of the flushed inlet actuator. As the water flow rate 

Q and S.E.S. velocity V came to a steady state, the inlet 

o 

throat area (A^) increased to its steady state value so that 

V^/V^ is exactly equal to 0.8. 

The pump torque ) and engine torque (^_) were 

p E 

plotted to check back flow of power during deceleration. 
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CONCLUSIONS 

Prom the result of the open loop simulation to a 
step fuel flow rate for the waterjet propulsion for S.E.S. 
several conclusions can be made. The first conclusion is 
that despite the pump being designed not to cavitate at 
hump speed with a thrust to resistance ratio of 1.2, cavi- 
tation can still occur during transients since designed 
pump head and designed pump flow rates for the hump speed 
can be reached way earlier before the S.E.S. goes to hump 
speed during acceleration. However, such cavitation only 
occurs for a few seconds. 

The second conclusion is that for a fixed designed 
area nozzle the pump is operating very close to the designed 
<p and ^ at steady state for all power levels; likewise, 
the gas turbine is operating very close to the point of 
maximum efficiency for all power levels at steady state. 

The third conclusion is that in the criteria for 
designing the speed at which the inlet area can change, there 
are two variables involved, namely the flow rate Q and the 
S.E.S. velocity . During transients it was noted that Q 
changes faster than Vq. Thus, if local cavitation is to 
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be avoided in the inlet during transient, the rate of 
change of the term Q/A^V^ should be considered. 

As already discussed in the MBODE analysis of poles 
and zeroes locations of the system the prime mover and the 
pump dynamics can be decoupled from the S.E.S. motion 
dynamics so that the design of a stable control for pump 
RPM can be simplified. A similar finding was found to be 
true also between the inlet area dynamics and the S.E.S. 
motion dynamics. 

In the fuel demand scheduled rate limiter incorporated 
in the RPM control of piamp, cavitation can be avoided during 
transient by setting the rate limiter to some lower value 
for the normal operation of the S.E.S. For the hard accel- 
eration maneuvers the fuel demand scheduled rate limiter 
can be set according to what the gas turbine characteristic 
may allow for surge and overtemperature . 

The control algorithm of Fig. 39 works for a rate 
limited and hard limited actuator and is recommendable 
since analog or digital implementation is relatively simple. 
It also works throughout the whole nonlinear operating range 
of the water jet system. 
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RECOMMENDATIONS 

The paper presented as it is, is like a textbook 
approach to the basic principles involved in the system. 
Several engineering approximations were made for a sim- 
plified digital simulation and control of the pump RPM and 
flushed inlet area. The following recommendations are 
then made for future work in this area of study. 

1. In the control of the flushed inlet, only the 
throat area was considered in this work. This was in- 
tended to serve as a reference in the actual reshaping of 
the whole flushed inlet geometry consistent with maintaining 
maximum possible pressure recovery, avoidance of cavitation 
and minimum drag. The mechanical implementation of this 
scheme deserves some further study. 

2. The fuel demand scheduled rate limiter still 
needs some refinement to avoid cavitation during acceleration 
for S.E.S. normal operation. Such a limiter should be 
consistent with the operating height of the pump above the 
waterline, and S.E.S. motion, prime mover and pump dynamics. 

3. In order for the S.E.S. to maintain a straight 
line motion there is a need for another set of controls to 
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regulate the thrusts of all the water jet systems. Similarly, 
turning maneuvers with and without the use of thrust de- 
flectors deserve some further study. 

4. Backing maneuvers with thrust reversers were 
not included in the study due to the absence of data on 
flushed inlet performance in the backward motion. This 
problem should also be investigated. 

5. It will be recommendable to implement a hybrid 
simulation of the system so that the complete gas turbine 
model as presented in Chapter II can also be included. 
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APPENDIX I 

The following are the particulars of the S.E.S. 
under study: 

Length to beam ratio =2.37 

Cushion pressure to length ratio =2.5 

, 2 

Cushion pressure = 412 Ibf/ft 

Displacement = 2109 tons 

Number of lift engines = 2 

Number of propulsion engines = 6 

Designed cruise speed = 80 Knots 

Maximum intermittent engine rating = 22,500 HP 

Maximum continuous engine rating = 20,500 HP 

Range = 2,000 miles 

Hump speed = 30 knots 

Bare hull resistance (smooth sea, cushion borne): 



Velocity 

(Knots) 


Resistance 

(Ibf) 


Velocity 

(Knots) 


Resistance 

(Ibf) 


28.45 


422918 


45 


299578 


30 


429438 


50 


282030 


35 


367391 


55 


271360 



40 



326606 



60 



265899 
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Velocity 


Resistance 


(Knots) 


(Ibf) 


65 


261984 


70 


258695 


75 


259040 


80 


263038 



Velocity 


Resistance 


(Knots) 


(Ibf) 


85 


269772 


90 


278874 


95 


290098 


100 


303252 
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APPENDIX II 



The steady operating point values of S.E.S. velocity 



(V^^) , pump 


RPM (N ) , 
o 


water flow rate per 


pump 


) , 

1 


throat area 


per inlet 


(A- ) and 
o 


reduction 


gear efficiency 


(W^ ) are; 
■‘■o 


different 


values of 


fuel flow 


rate per 


engine 


"fo 

(Ibm/hr) 


^oo 

(ft/ sec) 


^o 

(RPM) 


Qo 

(f t^/ sec) 


-^io 

(ft^) 


^ RGo 


800 


9.3 


109 


47,38 


6.307 


0.858 


2,000 


20.1 


241.8 


101.6 


6.318 


0.92 


3,000 


28.96 


346.8 


146.4 


6.319 


0.944 


4,000 


34.37 


410.3 


173.7 


6.317 


0.952 


5,000 


43.22 


465 .6 


197 .5 


5.317 


0.957 


6,000 


46.62 


514 


218.0 


6.317 


0.961 


7,000 


49.67 


553 .9 


235.6 


6.317 


0.962 


8,000 


103 .0 


589.9 


251.2 


6.317 


0.963 


8,500 


124.0 


605 .8 


263 .9 


3.202 


0.964 


9,000 


124.0 


620.9 


280.4 


2.822 


0.965 


9,500 


135.2 


635.6 


292.0 


2.701 


0.966 


10,000 


142 


647 .9 


298.6 


2.628 


0.966 
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From the above operating point values, power turbine 

torque ( 't" ), pump flow coefficient ( ) and head coef- 

Eq ° 

ficient { ) can be determined. Thus the matrix A of the 

' o 

linearized waterjet system equation (3-8) can be evaluated 
so that: 



for =800 

fd 



A = 



0.0264 


o 

• 

o 


0.00514 


o 

• 

o 


o 

• 

o 


o 

• 

o 

1 


0.0 


-0.0293 


-0.00757 


0.001302 


o 

• 

o 


o 

• 

o 


0.5172 


1.25 


-3.024 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


-4.459 


o 

• 

o 


-0.5 


4.0 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


O 

• 

o 


o 

• 

CN 

1 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


-0.5 



for = 5,000 



A = 



0.1088 


o 

• 

o 


0.02139 


O 

• 

o 


o 

• 

o 


1 

o 

• 

o 


o 

• 

o 


-0.0106 


-0.04479 


0 .001449 


o 

• 

o 


o 

• 

o 


2 .119 


5.2 


-12.337 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


-20.179 


o 

• 

o 


1 

o 

• 


3.4 


o 

• 

o 


0 .0 


0.0 


0.0 


0.0 


0 

• 

(N 

1 


o 

• 

o 


o 

• 

o 


o 

• 

o 


0.0 


o 

• 

o 


o 

• 

o 


-0.5 
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for W = 10,000 



- 0.03283 


0.0 


0.02556 


0.0 


0.0 


- 0.1242 


0.0 


- 0.1754 


0.0019 


0 .0014 


0.0 


0.0 


7 .7 


7 .704 


- 19.706 


o 

• 

o 


O 

• 

o 


0.0 


o 

• 

o 


-32 .506 


o 

• 

o 


- 0.5 


2.65 


O 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


O 

• 

o 


1 

to 

• 

O 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


o 

• 

o 


- 0.5 



The above constants were used in designing controller 
gains for small disturbance. 



I 

I 




118 



APPENDIX III 

Subroutine EQSIM is required by DYSYS for the solution 
of the dynamic equation of the system. The DYSYS program is 
available at M.I.T. mechanical and civil engineering computer 
facilities . 

Subroutine LM250 represents the LM2500 performance 
map of Fig. 8. The inputs to the subroutines are power 
turbine RPM and fuel flow rate and the output is power 
turbine torque . 
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